Instability of Non-vortex State toward a Quantized Vortex in Bose-Einstein 

Condensate under External Rotation 



Tomoya Isoshima* and Kazushige Machida 
Department of Physics, Okayama University, Okayama 700-8530, Japan 

(February 1, 2008) 



OS 

(N 



o 

CO 

-4-* 



o 
o 



> 

(N 
O 

o 
m 
o 

On 
o3 



i 

C 

o 
o 



X 
S3 



The instability condition of the non-vortex state toward 
vortex formation is examined within the Bogoliubov theory 
when a Bose-Einstein condensate is under externally forced 
rotation. The obtained critical angular velocity combined 
with the previous stability conditions for a votex yields a de- 
tailed phase diagram in the critical velocity vs the system 
parameter. This facilitates vortex formation experiments for 
alkali atom gases confined in a harmonic potential. 

PACS numbers: 03.75.Fi, 67.40.Vs, 67.40. Db, 05.30.Jp. 



The quantized vortex is a manifestation of macro- 
scopic quantum phenomena in quantum fluids and have 
been observed in superconductors and superfluid helium. 
Since the experimental realization jjj of Bose-Einstein 
Condensation (BEC) in alkali-atom gases, the vortex 
problem in weakly interacting Bose gas has been inves- 
tigated theoretically M. So far there is no experimental 
report to find a vortex in the present BEC systems and 
several trials are ongoing. 

It is well known that the Gross-Pitaevskii (GP) theory 
P,|l is valid for describing a weakly interacting Bose sys- 
tem which is the case for the present alkali-atom gases of 
Rb and Na. Based on the GP theory the single vortex 
state is discussed by several authors j4|j| . More recently 
Butts and Rokhsar || describe a many vortex state in 
rotating BEC, giving rise to an interesting vortex config- 
uration with many vortices. This may be directly observ- 
able experimentally. The time evolution of vortex is also 
studied by a time dependent form of the GP equation . 

To investigate further fundamental issues on a quan- 
tized vortex, we resort to the Bogoliubov theory 
This microscopic framework can describe the conden- 
sate and the non-condensate components simultaneously, 
yielding the collective excitation spectrum by solving an 
eigenvalue equation ]lQ - |l5|l. 

In our series of papers |lfj-jl2| , we have been studying 
the stability problem of an isolated singly quantized vor- 
tex. The main conclusions drawn from our studies are: 
(1) Contrary to the conclusion in classical papers ||f§, 
BEC cannot sustain the stable vortex in an infinite sys- 
tem at T=0 |T(| . This conclusion coincides with Rokhsar 
(2) The finite temperature which increases the non- 
condensate component relative to the condensate one ef- 
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fectively stabilizes it in a harmonically confined system 
p| . (3) A pining potential added to a harmonic con- 
fined potential also stabilizes the quantized vortex 
We note that concerning the conclusion (2), Stringari flq | 
gives the phase diagram of the stability in a wide range 
of T by comparing total free energies. 

In |l2] ], apart from the above intrinsic vortex stabil- 
ity problem, we have discussed the vortex stability under 
externally forced rotation at T = 0. There we introduce 
two different criteria for stability: The global stability is 
defined by comparing the two total energies of the non- 
vortex state and the vortex state. The other local stabil- 
ity is defined by the intrinsic stability condition for the 
vortex itself, evidenced by the positiveness of the lowest 
eigenvalue. According to these criteria, we obtain the 
phase diagram (see Fig. 7 in Ref. |12|) in the rotation an- 
gular frequency Q vs the system parameter an z (a: the 
scattering length and n z : the particle density per unit 
length), which shows the three divided regions upon in- 
creasing il: (A) the vortex unstable region, (B) the single 
vortex stable region, and (C) the single vortex unstable 
region where many vortices may appear. 

In order to quantify the problem further and refine this 
phase diagram, we investigate the local stability of the 
non-vortex state under external rotation. Through this 
study we may gain a valid picture of the actual nucleation 
process when we begin to rotate a BEC quite slowly from 
zero, which has not been examined before. It is our hope 
that this study fosters the observation of a vortex in the 
present BEC systems. 

Since the formulation and computational procedures 
are same as before [l^] , we briefly introduce the main 
equations and the notations. Let us start the fol- 
lowing Hamiltonian for interacting Bosons with the 
s-wave scattering length a: H rot = H — oj ■ J r x 
p(r)dr where H = J dr*t(r){-^^ + V(r) - ^}*(r) + 
^dr^t(r)^(r)^(r)^(r), where p(r) is the momentum 
operator, /i is the chemical potential, and V(v) is the 
confining potential. The mass of a particle is m. 

Following the standard method, we decompose the 
field operator * as *(r) = ip(r) + 0(r), 4>(r) = (*(r)). 
We substitute the above decomposition into the above 
and ignore the higher order terms such as ip^ipip, ijfi <\\) , 
and ip' ip' ipip terms. Then we rewrite the operator ip( r ) as 
= ^ q [u q (r)r] q — v*(r)?yj] where rj q are the annihi- 
lation operators, u q (r) and v q (r) are the wave functions, 
and the subscript q denotes the quantum number. 

We consider a cylindrically symmetric system and a 
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vortex line, if exists, passes through the center of a cylin- 
der, coinciding with the rotation axis. We use the cylin- 
drical coordinates: r = (r, 9, z). The system is trapped 
radially by a harmonic potential V(r) — im(27ri/) 2 r 2 and 
is periodic along the z-axis whose length is L. 

We focus on the lowest eigenstates of the momen- 
tum along z-axis and the quantum number q is writ- 
ten as (q r ,qe), where q r = 1,2, ••• and qg = 0, ±1, 
The condensate wave function (f>(r) is expressed as 

cf>(r,9,z) — cf)(r)e lwS . (0(r):real) In similar way, u q (r) = 
Uq ( r y{ qe +w)e v ^ T j _ Vq (r) e ^ie-w)e _ The cage w = q 

corresponds to a system without a vortex and w — l(w = 
2) with a singly (doubly) quantized vortex. 

The condensate has N particles in the length L and 
the area density is n z = N/L. We introduce the den- 
sity unit rio = yjhv / g where g = -iTrh 2 a/m, the length 
unit \Jlmhv /%. The energy (e', //) are scaled by hv. 
The wavefunctions (</>, u, v) are scaled by y/nQ. The an- 
gular velocity w is also scaled by 2i:v and the normalized 
rotation is £1 = tt—. We use the normalized values of 
fi, s, <j), u, v, r, and fl below. 

The condition that the first order term in ijj(r) of our 
Hamiltonian vanish is the Gross-Pitaevskii equation 

d?_ ld L _ vfi\ _ 
dr 2 r dr r 2 J ^ 
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The normalize condition is J \<fi(r)\ 2 r dr 
condition that the Hamiltonian be diagonalized gives the 
coupled eigenvalue equations for u q (r), and v q (r) whose 
eigenvalue is e q : 
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+ — + 2<p 2 (r) + Q(w + qe)\ u q (r) 

-4> 2 {r)v q {r) = e q u q (r), 
d 2 1 d (w — qe) 2 
dr 2 r dr r 2 
r 2 

+ — + 24> 2 {r) + VL{w - q e )\ v q (r) 

~(f) 2 {r)u q (r) = -e q v q (r). (3) 

The normalization conditions are 

/ { u p( r ) u q( r ) - v p( r ) v q( r )} rdr = Y;5 P , q - We determine 
the signs of qg and e by the above normalization condi- 
tion. The set of the Eqs. (|l|), (@), and (||) constitute the 
Bogoliubov theory of our system. 

The contributions of the non-condensate to the total 
energy is negligibly small and we will use the contribution 
from the condensate 
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l£ - *rW) + T\Hr)\ 2 + W)\ A ]rdr - wQ as the total 
energy of the system. Here 4an z is a dimensionless num- 
ber. Note that the coefficient of u> in Ed, is a constant. 



Since the parameter ^ in the conditions of normalization 
does not change the following results within the present 
Bogoliubov framework, the system is characterized by 
the number an z and the rotation f2. 

We solve the coupled equations; Eqs. (Q), (||), and (||) 
for the gas of 23 Na atoms trapped radially by a harmonic 
potential of j/Hz. The area density per unit length along 
the z-axis is chosen to be an z = 2.75 ~ 137.5. The 
density profile of the condensate cj) 2 (r) is solved from the 
GP equation Eq. ([!]) . Figure [j] shows the profiles of the 
system with w = 0, 1, and 2. 
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FIG. 1. The particle number density distribution <^> 2 (r) 
along the radial direction. Dashed (dotted) line is the singly 
(doubly) quantized vortex case w = 1 (w = 2) while solid line 
corresponds to the non-vortex case. To draw this figure in 
actual units, we use the following parameters: the scattering 
length a — 2.75nm, m = 3.81 x 10~ 26 kg, the radial trapping 
frequency v = 100Hz, and n z = 2 x 

The eigenvalues e's evaluated from Eqs. (|^) and (j^) 
should be all positive relative to the condensate energy 
at zero for the system to be stable. This defines the 
local stability of the system. In fl2|| , we have derived 
the critical fi's of the local stability in vortex state. The 
lower criterion is £^ ca i and the upper criterion is f^ cal - 
For example, Qf ocal = 0.16085 and fjU cal = 0.44610 when 
an z — 55. 

Also the non- vortex state, which has the winding num- 
ber w = 0, could become locally unstable under the large 
external rotation. In other words, some excitation level 
Eq at larger qg can become negative. These negative ex- 
citations correspond to "surface excitations" derived by 
Dalfovo @. 

Here, we calculate the critical Q in the non- vortex state 
to refine the phase diagram (Fig. 7 in @) of stabili- 
ties. In Fig. H the lowest edge of the eigenvalues in the 
non- vortex state is displayed for selected f2's. There is 
a critical value of £l w =o above which some eigenvalue(s) 
in qg > become(s) negative, signalling an instability of 
the non-vortex state. When il < fl w —o all the eigenval- 
ues are positive. When an z — 55, ^^=0 = 0.41145. The 
resulting curve for tt w= o is drawn as a function of an z in 
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FIG. 2. The lowest edge of the eigenvalues along qg for 
selected angular velocities Q's. Circles, filled circles, trian- 
gles, and filled triangles correspond to Q = 0, 0.15, 0.3, and 
0.45 respectively (an z = 55). Another eigenvalues is densely 
distributed at each qg. The eigenvalues in qg < move up 
while the eigenvalues in qg > move down as Q increases. 
At larger SI some positive eigenvalue of qg around qg ~ 10 
becomes negative as fi increases. The eigenvalue at qg = 
does not change with Q. 
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FIG. 3. The four main critical velocities f2's; Sl local , 
^global; Qw=o, and f2}o Ca i are plotted in solid lines as func- 
tions of an z . These four lines divide the whole area into five 
regions (i) - (v). Each of regions are explained in text and in 
Fig. |^. The dashed line denotes f2g[ obal . The normalization 
condition for Eq. (^); J |0(r)| 2 rdr = 4an z will be helpful to 
understand the value an z . 

Besides the above local stability conditions for the non- 
vortex and vortex states, there is the global stability con- 
dition for fi. The comparison of the two energies in 
the systems with and without a vortex defines the "global 
stability" . The critical value between non- vortex state 
and single vortex state is Slg loba i an( i that between single 
vortex state (w = 1) and double vortex state (w = 2) 
is ^global- The critical value is Slg lob ai = 0.22215 and 
^global = 0.42071 when an z = 55. 



As mentioned previously, the system is characterized 
by the number an z and the rotation SI. In Fig. || we 
summarize the four critical Si's; 0}^^, Slg lobal , ^w=o, 
and fi£ ca j as a function of an z . We can divide the whole 
range of SI into five regions, which are labeled as (i) to 
(v) in Fig. |. Each region is characterized in Fig. ^ where 
the schematic energy curve for each region is shown. 

As SI increases from zero, 

(i) the non-vortex state is stable both locally and glob- 
ally for < Q < fifc^. 

(ii) The singly quantized vortex also becomes stable lo- 
cally for SlJ ocal < SI < Sl^j obal . But the vortex state has 
higher energy than that in the non-vortex state. There- 
fore, a vortex does not nucleate spontaneously in this 
region. 

(iii) The energy of this vortex state becomes lower than 
that in the corresponding non-vortex state for Slj; lobal < 
SI < Q w —o. As it is seen from Fig. ^, the vortex state 
is locally stable now, but the non-vortex state is still lo- 
cally stable. Therefore, the non-vortex state must tunnel 
quantum-mcchanically through an energy barrier to the 
vortex state. 

(iv) The vortex state can nucleate spontaneously in 
£l w =o < < ^locai' because the non- vortex state is un- 
stable both locally and globally, and there is no reason 
for the system to stay in the non-vortex state anymore. 
A vortex is now bound to appear. 

(v) Since single vortex state becomes unstable, many vor- 
tices may appear in this region: Sl^ ca] < S7. 

It is to be noted that even above the £lgi b a i which is 
usually regarded as a critical angular velocity of vortex 
formation, the non-vortex state remains locally stable in 
region (iii). The upper edge of (iii) region (fiu, = o) ex- 
ceeds twice of ^gi bai a t larger an z . This means that the 
angular velocity much larger than Sl^ lobal is needed to 
create a vortex. 

The comparison between E$^ w= i and E^^—i defines 
^global- The energy E^w—2 may become lower than that 
of the state with two singly quantized vortices (this en- 
ergy cannot be calculated in this axial symmetry formu- 
lation). But in the discussion below, we treat S^gi b a i as 
the (tentative) angular velocity to stabilize many vortices 
states. 

We have ignored Slg lobal in Figs. || and f|. This does 
not contradict with the above explanations in (iii) or (iv) 
region. But at much smaller an z , the ^gi b a i appears 
in (v) region. It means that when the w = 1 state is 
unstable locally toward many vortices state, the energy of 
w = 2 vortex state is higher than that of the w = 1 state. 
It seems to contradict with the illustration in the first line 
of Fig. ||. A possible explanation is as follows: The radius 
of w = 2 vortex is so large that it is comparable with the 
system's radius with w = 1. This increases E w= 2 and 
thus SlU Qbal . 
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FIG. 4. From left to right column, the ranges of Q, the 
schematic energy configuration near non-vortex and vortex 
states, the local stability, and the global stability. The number 
of the last column corresponds to that of Fig. pi 

In order to obtain more insights into the vortex stabil- 
ity problem in the Bose-Einstein condensation of alkali 
atom gases confined in a harmonic potential, we have ex- 
tended our previous work p||ll| to the case where the 
cylindrical system is under forced rotation. 

Our calculation is done within the framework of the 
non-selfconsistent Bogoliubov theory and yielded the vor- 
tex stability phase diagram shown in Fig. ||. The stabil- 
ity of the vortex state is examined by the two different 
ways. They are the local stability and the global stability. 
The local stability is evaluated both in single vortex state 
and in non-vortex state. These stability conditions yield 
the whole phase diagram of the vortex stability problem 
shown in Fig. || In addition to the usual critical angu- 
lar velocity figiobaij we quantitatively obtain the values 
of £l w =o and g i o b a i which describes the instability of the 
non- vortex state toward the vortex state. 

The previous cases [jlO|,[ll]], which are consistent with 
the present results, correspond to the £1 = axis in Fig. || 
where the vortex is intrinsically unstable. This instability 
is seen to occupy a finite range of Q, not only the isolated 
line confined at fl = 0, for various densities n z and the 
scattering lengths a as another axis. It is found that the 
angular velocity much larger than ^ iohai is needed to 
create a vortex. With the previous finite temperature 
calculations [ pT| , the present calculation concludes that 
alkali atom Bose gases in BEC can sustain and exhibit 
the stable vortex in certain range of temperature and 
external rotation. 



There are many problems left concerning the vortices 
in BEC. About the appropriate value of angular velocity, 
we need the calculations in full three dimension geome- 
try. The process to gain/lost the vortex core will be de- 
scribed better if we can treat the time evolution of both 
the condensate and the non-condensate . The way to 
give the angular momentum to the condensate (or cloud 
of atoms before condensation) is the important problem 
for experimentalists p!9[ . It is our hope that this study 
fosters the observation of a vortex in the present BEC 
systems. 

The authors thank L. Pitacvskii for raising an adequate 
and timely question regarding to the instability of the 
non-vortex state, which leads to the present work. 



[2 

[3. 

[4] 
[5] 
[6] 
[7] 

[8 
[9 

[10 

[11 

[12 

[13 
[14 

[is; 



See for review: F. Dalfovo et al 



e-prmt 
Mod. Phys 



cond- 



71, 



mat/9806038| [To be published in Rev. 

20, 454 (1961); J. Math. 



April (1999)]. 

E. P. Gross, Nuovo Cimento 
Phys. 4, 195 (1963). 

L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 646 (1961) [En- 
glish Transl. Sov. Phys. -JETP 13, 451 (1961)]. 
M. Edwards et al, Phys. Rev. A 53, R1950 (1996). 
S. Sinha, Phys. Rev. A 55, 4325 (1997). 
D. A. Butts and D. S. Rokhsar, Nature 397, 327 (1999). 
B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. 
Lett. 80, 3903 (1998). 

N. Bogoliubov, J. Phys. (Moscow) 11, 23 (1947). 
A. L. Fetter, Phys. Rev. 138, A709 (1965); 140, A452 
(1965); Ann. Phys. (N.Y.) 70, 67 (1972). 
T. Isoshima and K. Machida, J. Phys. Soc. Jpn. 66, 3502 
(1997). 



T. Isoshima and K. Machida, e- print cond-mat/9807250 
[to be published in Phys. Rev. A 59, 1 March (1999) 



T. Isoshima and K. Machida, e-print cond-mat/9811367 
[to be published in J. Phys. Soc. Jpn. 68, No.2 (1999)] 
R. J. Dodd et al., Phys. Rev. A 56, 587 (1997). 
D. S. Rokhsar, Phy s. Rev. Lett. 79, 2164 (1997); e-print 



cond-mat/9709212 



A. A. Svidzinsk y and A. L. Fetter 
(1998); e-print |cond-mat/981134£ 



Phys. Rev. A 58, 3168 
A. L. Fetter, e-print 
cond-mat/9808070| [To be published in J. Low Temp. 

Phys-]- 

S. Stringari, e-print |cond-mat/9812362 . 

F. Dalfovo et al, Phys. Rev. A 56, 3840 (1997). 

T. Nikuni, E . Zaremba, and A. Griffin, e-print cond 



mat/9812320 



[19] K.-P. Marzlin and W. Zhang, Phys. Rev. A 57, 3801 
(1998); 57, 4761 (1998). 



4 



